Experimental Astronomy manuscript No. 

(will be inserted by the editor) 



Simulations of Galactic Cosmic Ray Impacts on the 
Herschel/PACS bolometer Arrays with Geant4 Code 

C.Bongardo ■ P.M. Andreani ■ K.Okumura • 
B.Horeau • A. Claret • G.De Zotti ■ 
R.Giannitrapani • F.Longo 



Received: xx/yy/200z / Accepted: xx/yy/200z 



Abstract The effects of the in-flight behaviour of the bolometer arrays of the Her- 
schel/PACS instrument under impacts of Galactic cosmic rays are explored. This in- 
strument is part of the ESA-Herschel payload, which will be launched at the end of 
2008 and will operate at the Lagrangian L2 point of the Sun-Earth system. We find 
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that the components external to the detectors (the spacecraft, the cryostat, the PACS 
box, collectively referred to as the 'shield') are the major source of secondary events 
affecting the detector behaviour. The impacts deposit energy on the bolometer chips 
and influence the behaviour of nearby pixels. 25% of hits affect the adjacent pixels. 
The energy deposited raises the bolometer temperature by a factor ranging from 1 to 6 
percent of the nominal value. We discuss the effects on the observations and compare 
simulations with laboratory tests. 

Keywords Instrumentation: detectors, bolometers • Galaxy: cosmic rays • ISM: 
cosmic rays 

PACS PACS: 95.55.Rg ■ PACS 96.40.De 
1 Introduction 

The European Space Agency Herschel satellite will be launched at the end of 2008 with 
an Ariane-5 rocket and will operate at the Lagrangian L2 point of the Sun-Earth system 
(see ESA web page: www.rssd.esa.int/SA-general/Projects/Herschel). Herschel is the 
ESA fourth cornerstone mission and will perform imaging photometry and spectroscopy 
in the far infrared and submillimetre part of the spectrum. The Herschel payload 
consists of two cameras/medium resolution spectrometers (PACS and SPIRE) and a 
very high resolution heterodyne spectrometer (HIFI). 

The PACS (Photo-conductor Array Camera and Spectrometer) instrument will 
perform efficient imaging and photometry in three wavelength bands in the range 60- 
210/im, and spectroscopy and spectroscopic mapping with spectral resolution between 
1000 and 2000 over the same wavelength range or short segments. 

PACS is made of four sets of detectors: two Ge:Ga photoconductor arrays for the 
spectrometer part and two Si-bolometers for the photometer part. On each instrument 
side each detector covers roug hly half of the PACS bandwidtlfl 

It is well known that cosmic rays may influence strongly the detector behaviour in 
space. The performances of Infrared Space Observatory (ISO) detectors were largely 
affected for time periods long enough to corrupt a large amount of data [T]. Our goal 
is to exploit our present knowledge about the detectors and the cosmic environment to 
understand how the detector behaviour changes and how we may retrieve the lost infor- 
mation and/or extract the astronomical signal from the sources through an appropriate 
data analysis tool. 

In this paper we focus our attention on the Si-bolometer arrays of the PACS in- 
strument. These detectors differ from those on board of ISO as they are not of photo- 
conducting type and very likely less affected by cosmic ray hitting. In a previous paper 
we dealt with the Ge:Ga photoconductors of the same instrument ([2]). 

The paper is organized as follows: in Sections [2] through [5] we briefly describe the 
used simulation toolkit, the input detector design, the Galactic cosmic ray spectra and 
the physics list. In section[6]we discuss very simple tests to check the reliability of the 
present simulations. The results are reported in section [7] and discussed in section [8] 
while in section we compare simulations with laboratory tests. Conclusions about the 
effect on observations are to be found in 1101 

More about PACS can be found in the PACS Web pages: pacs.mpe-garching.mpg.de and 
pacs.ster.kuleuven.ac.be 



Fig. 1 Left: a highly magnified picture of one pixel (the Si grid) of size 0.75mm. Such pixels 
are assembled into 16x16 pixel subarrays (middle), which in turn are assembled into the focal 
plane matrices. Right: a picture of the focal plane PACS blue bolometer camera. 

2 The Geant4 Monte Carlo code 

Monte Carlo simulations were carried out with the Geant^EI (hereafter G4) toolkit, 
which simulates the passage of particles through matter. G4 provides a complete set 
of tools for all the domains of detector simulation: Geometry, Tracking, Detector Re- 
sponse, Run, Event and Track management, Visualization and User Interface. An abun- 
dant set of physics processes handle the diverse interactions of particles with matter 
across a wide energy range, as required by the G4 multi-disciplinary nature; for many 
physics processes a choice of different models is available. For any further information 
see the Geant4 Home page at: geant4.web.cern.ch/geant4/. 

3 The Detector Design 

The PACS camera is made of two bolometers: one called 'blue' working at short wave- 
lengths (60-130/im) and one 'red' at long wavelengths (140-200/im). The blue one con- 
sists of 8 chips (4x2), each of dimensions of 15.78 x 12.63 mm and comprising 16 x 16 
pixels. The chips are separated from each other by 120 /im and have a projected field 
of view in the sky of 1.75' x 3. 5'. The red one is made of 2 chips covering the same sky 
area, each of dimensions of 12.63x12.63 mm and comprising 16 x 16 pixels. 

The detector module of the PACS bolometers is made of different parts: each pixel 
consists of an absorber supported by a silicon grid 5/im thick with a thermometer, 
and is suspended in an Inter-Pixel- Wall (IPW) grid 450/im thick which has a reference 
thermometer and acts as a thermostat, of aluminum, Titanium and Silicon layers, 
of a Silicon strip, and of reflecting slabs [3|. The output signal is a function of the 
temperature difference between the pixels and the reference thermometers. 

Figure Q] shows a picture of the real PACS blue bolometer detectors, while Figure 
[2] shows how we have modeled their geometry. This last configuration is that used for 
for the simulations with the G4 software. Details of the modeled geometry are given 
in Appendix A. 

The ensemble of all the components 'external' to the detector, i.e. the PACS box, 
the cryostat and the spacecraft is hereafter called 'shield'. The dimensions of this 

2 Every simulation was run with the 7.0 patch 01 version, vith CLHEP 1.8.1.0 with a gcc 
3.2.3 compiler 



Fig. 2 The PACS blue bolometer array. Grey: the Al layers. Green: the Ti layers. Blue: the 
Si layers. Cyan: the Si strip. Black: the IPW. 



latter have been chosen to be big enough to contain the PACS detectors. A detailed 
simulation should in principle include a model of the exact geometry of the spacecraft. 
For the PACS bolometer, however, we approximate the spacecraft, the cryostat and 
the telescope with an aluminum sphere, 11 mm thick and at a temperature of 80 K, 
which is the nominal temperature of the passively cooled Herschel telescope. The inner 
radius of the 'shield' is 20 cm and the outer radius 21.1 cm. This aluminum sphere 
corresponds to 15.7 Kg, i.e. 5 percent of the telescope mass. We have carried out the 
same tests described in this paper by varying the telescope temperature from the 80 
K value and it turns out that the results do not change. 

The G4 software then requires the design of a world volume containing all the 
detector parts and the specification of the environment in which the detector is placed. 
We made the world volume a little bit larger than the 'shield', i.e. a sphere with a 21.2 
cm radius. 
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Table 1 Rates (number of nuclei per second hitting a sphere of 21.2 cm radius) for different 
GCR types and solar activity levels. Values reported here are the integrals of Figure 3. 



GCR Type 


Solar min 

rates (#/s) 


Solar max 
rates (#/s) 


H 


6639 


2452 


He 


649 


270 


Li 


4.4 


1.7 


C 


19.7 


8.2 


N 


5.3 


2.2 


O 


18.4 


7.7 



4 Galactic Cosmic Rays spectra 

Herschel will be injected in a orbit around L2. The Galactic Cosmic Ray (hereafter 
GCR) rate in the L2 environment has been measured by the WMAP satellite [3] and 
found to be very close to that of the geostationary orbit, for which the GCR spectra 
are known. We use the Cosmic Ray Effects on Micro-Electronics (CREME) model set 
both for the solar-minimum (GCRs maximum) and for the solar-maximum environ- 
ment [5]. CREME treats geosynchronous and near-Earth interplanetary environments 
as identical, there is no geomagnetic shielding and no trapped particles Q The CREME 
interplanetary flux model (galactic cosmic rays, anomalous cosmic rays, and solar en- 
ergetic particles) is based on measurements at Earth (1 AU). CREME uses numerical 
models for the ionizing radiation environment in near-Earth orbits to evaluate the 
resulting radiation effects on electronic systems in spac^3. 

Our simulations take into account the most common GCR particles in the near- 
Earth heliosphere: protons, alpha particles and nuclei of Li, C, N, O (see Figure[5J. We 
had to neglect heavier ions as physical processes in Geant4 lack a detailed comparison 
with experimental data and reliability. Fe ion, although with an abundance lower of 
two order of magnitudes than that of Oxygen may cause a big jump in the detector 
temperature because of its larger mass. The current version of G4 is not reliable to 
follow these processes and they are then neglected in this work. In Table [1] we report 
the rates (number of particles per second) of each GCR type considered, computed by 
integrating the fluxes shown in Figure [3] over the 'shield' surface and particle energy. 
We have considered here that the detector is hit by the flux coming from all directions 
and that a Herschel will reach the L2 point at a solar minimum but it will operate until 
an epoch of increasing solar activity. Both phases of solar activities are thus considered. 



4.1 General Particle Source 

The tool is built in such a way that once geometry and physics models are chosen and 
set up it can be driven via a scripting User Interface (i.e. a text macro file) in which 
the user specifies: (a) the graphical visualization and/or (b) the type and the number 
and energy of the incident particles (only part of which will reach the detector). One 
single particle with a determined energy represents an event. 

3 Strictly speaking, this is not correct, but the differences arise at energies below ~ 5 
MeV/nuc, which are irrelevant for single-event effects. 

4 In this work we used the CREME86 version, which was more extensively tested 
and compared with measurements. An updated version, CREME96, is available at: 
https:/ / creme96. nrl. navy.mil/ 
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Fig. 3 Differential GCR Spectral Energy Distributions (SEDs) at the solar minimum (left) 
and maximum (right). Each curve corresponds to a type of particle; from top to bottom: 
protons, a, C, O, Li, N. Values are from CREME86. 

We simulated the GCRs seeding via the General Particle Source (GPS; [B], [7], 
[S]), a code for the spacecraft radiation shielding, that is part of the G4 radiation 
transport code. The user must specify the input parameters of the particle source: (1) 
the spatial and (2) the energy distribution. Every time G4 requires a new event, the 
GPS will randomly generate it according to the specified distributions. We adopted an 
angular distribution cosine law on a spherical surface (the vacuum sphere), without 
any preferred particle momentum. The energy distribution follows the SEDs shown in 
Figure [3] (see also section [6j . 



5 The Physics Model 

We also had to define a list of physical interactions with the 'shield' and the different 
components of the bolometers. The model called by our physics list depends on the 
energy of the incoming particles. The reader may find a detailed list of the physical 
processes together with our own choices in [2] and in Appendix B. 



5.1 Thresholds 

In order to avoid the infrared divergence some electromagnetic processes require a 
threshold below which no secondary particles are generated. Each particle has an asso- 
ciated production threshold, i.e. either a distance, or a range cut-off, which is converted 
to an energy for each material, which needs to be defined as an external requirement. A 
process can produce the secondaries down to the selected threshold (GEANT4 Users's 
Guide section 2.4.2). 

We have set different production thresholds (hereafter cuts) for each geometrical 
region. The G4 starting values, 0.7mm, could force a too large production of secondary 
particles for thick volumes such as the 'shield', and a too small production for thin 
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Table 2 The five regions of the bolometer camera and their specific cut values. 



Region 



Cut Value 



Thickness 



Shield 

Containers 

Si.Thick 

Si-Thin 

bolometer 



2 mm 
0.5 mm 

80 fim 
15 fim 

3 fim 



11 mm 

1.56 - 2.5 mm 
350 - 450 /im 
70 (im 
0.3 - 5 /im 



regions such as those of the IPW and the grid. We first ran our G4 code with the starting 
cut value of 0.7 mm for all volumes and then tuned them according to the bolometer 
geometry. We decided: 1) not to use very small cut values (below 1-2 ^m), since this 
could affect the validity of physical models at such small steps (fvantchenko, private 
communication); 2) to set cut values « 1/5 of the volume thickness. We identified 
groups of volumes with the same thicknesses: in particular we defined 5 different regions 
(see Table [2}. Starting from these assumptions we run simulations with five different 
cut values per volume, each one varying according to the thickness of the region under 
consideration. 

We find that the results do not depend strongly on the chosen cuts and the energy 
depositions for the different simulations were always consistent with each other to 
within la. 



5.2 Theoretical computation of the GCR flux 

fn order to disentangle the contributions from the different simulated parts of the 
instrument (detectors and 'shield') we evaluated the impact of the GCRs flux on the 
'naked' detector (without spacecraft shielding). In Table [3] we report the results from 
the simulated GCR flux hitting the detector only (results have been normalized to the 
frame rate fixed by electronic and software, 1/40 s). 

We compare them with the theoretical values computed as follows. Let $ be the 
particle flux (in cm -2 s _1 sr _1 ) and the unit surface dS, the number of impacts, n is 
given by (see i.e. [9]): 



where S is the total superficial area of the volume considered. This value (in Table|3]the 
'Th. Hit' column) should be considered as an upper limit, because it is the maximum 
theoretical value of the impacts number. 

These values (Okumura, 2004, internal report) are then compared with the average 
hit values obtained from the G4 simulations for the G4 starting cut values and for our 
cuts (see Table [3} . We did not take into account for this comparison the production 
of secondary events. Although statistics is poor, the results are in good agreement 
with the theoretical ones, with a slight better agreement in the latter case where the 
dispersion around the mean decreases. Table explains the kind of systematic uncer- 
tainty introduced by the calculation of the number of impacts per second as a function 
of the chosen cuts ("starting cuts" vs "chosen cuts"), i.e. the integration step size in 
comparison with geometric part size. 




sin 9 cos 9d8 = n$S, 



(1) 
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Table 3 Comparison of the theoretical estimate of the average hit rates (Th. Hit) with the 
G4 simulations with 'default' cuts and with our choice for them. Calculations were done for 
the blue bolometer without 'shield' and for two volumes and solar minimum values. 







Starting Cuts 


Chosen Cuts 


Volume 


Th. Hit (/s) 


{^impacts} / $ 


{^impacts} / 5 


Grid 


0.41 


0.5±0.3 


0.4±0.1 


IPW 


2.21 


2.5±0.9 


2.1±0.7 



Table 4 The n Bec /n pr i m ratio from simulations of GCR impacts on the bolometers. First 
column: GCR type; columns #2 and #3: average number of secondary events generated by 
a single primary event for the blue bolometer in the cases of solar min and max irradiation; 
columns #4 and #5: same as in columns #2 and #3, but for the red bolometer. 





Blue 




Red 






Solar min 


Solar max 


Solar min 


Solar max 


GCR 


(f^sec) /'Tlprim ('■ 


^sec) / Tlprim 




^sec) / Uprim 


H 


2.7±0.1 


3.6±0.1 


2.6±0.1 


3.5±0.3 


He 


6.6±0.2 


8.3±0.1 


6.6±0.2 


8.1±0.4 


Li 


11.1±0.4 


13.5±0.4 


11.0±0.2 


12.7±0.5 


C 


25.3±0.4 


33.6±0.7 


24.6±0.2 


33.0±0.4 


N 


32.3±0.5 


42.1±0.9 


31.7±0.7 


42.1±0.4 


O 


38.2±1.1 


51.6±0.5 


38.2±1.1 


51.5±0.5 



6 Preliminary Tests 

6.1 The n S ec/n pr i m ratio and production of secondary events 

GCR protons yield a broad variety of secondary products, from photons and electrons 
to unstable particles (i.e. 7r + ), to some heavier nuclei (i.e. 27 A1). For the other light 
nuclei, secondary events are mainly photons or electrons. Ion inelastic scattering may 
also generate pions, protons and electrons and lighter nuclei, but these events are much 
rarer. 

Secondary events are of utmost importance and are mainly produced in the 'shield'. 
Their production and their hitting on the detector chips make them the primary source 
of 'glitches' and largely affect the detector performance. 

Secondary events are originated also by the detector chips. As expected the blue chips 
produce more secondary events, since they are more massive than the red ones. But the 
difference between the two arrays is within the error-bars and then the main secondary 
events producer is still the 'shield'. 

TableUreports the number of secondary particles produced on average per primary 
event. We have listed all the secondaries produced without any energy threshold and 
independently of their energy release. 

The Tisec I '^prim ratio is higher in case of solar max, since, even if GCR.S are fewer, 
they are more energetic, that is they produce more secondary events. 
The risec/np r i m ratio is clearly increasing with the atomic number of the primary 
particle and the blue and red bolometers have same dependences, as the two cameras 
have the same thickness. Such a ratio is strictly dependent on the physics list we choose 
and, basically, on the goodness of the G4 models. G4 has been tested for C and lighter 
nuclei; therefore the ratios for N and O must be then taken with caution. The minimum 
proton energy required to produce secondary events is HOeV. 
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7 Results 

The bolometer arrays operate at 300 mK, the electronic readout system has a char- 
acteristic frame rate of 1/40 s. According to Table [1] we expect, in this time interval, 
about 166 proton impacts for a Solar minimum rate and about 61 at Solar maximum, 
and very low numbers for heavier ions. To get statistically significant results we should 
then carry out thousands of simulations, which would mean a too long computing time. 
We decided then to increase the number of particles and then normalize the outputs. 
We ran, then, five simulations for each nucleus, each one with 2xl0 5 particles. Results 
were normalized to the frame rate of the bolometer, dictated by software and electron- 
ics. This choice is further validated (see below) by the fact that the thermal effects are 
not cumulative (see below). 

7.1 Deposited energy 

We are particularly interested in the amount of energy deposited on the detector and in 
the number of pixels affected by one single GCR. The deposited energy is transformed 
into heat. The volumetric heat capacity of Si is cgj = 5.8 x 10~ 7 T 3 J/K/cm 3 . At 0.3 
K it becomes c Si = 1.57 x 10~ 8 J/K/cm 3 . 

The enhancement of temperature in the material is then related to the amount of 
deposited energy as follows: 

AE(J) = V ■ cgj • AT(K) = C Si x AT(K) (2) 

This equation implicitly implies that the efficiency of energy transfer is 100%. This 
is an ideal case and formally the resulting values must be considered as upper limits. 
However in view of the high energy involved in the process these values would not be 
too far from the real ones. As thermal effects are not cumulative in G4 we have added 
the deposited energies and then transformed them into temperature values. 
Impacts on the grid and on the IPW are identified and used to compute the deposited 
energy that will enhance the temperature of the bolometers, initially cooled to 300 
mK. Tables [5] and [6] report the mean enhanced temperature values (in mK) for both 
(blue and red) detectors for shielded simulations. 

7.2 Blue bolometer 

We consider first the effects on the blue bolometer array of the primary particle impacts 
on the 'shield' and the detector. The results are shown in Table [S] The number of hits 
are given per frame rate of the bolometer. 

A glitch event occurs when the energy deposited by the GCR exceeds some predefined 
energy threshold, which has not yet been defined for operation of the bolometers. We 
do not apply here any such criterion as an energy threshold cannot be defined yet. 
The numbers reported in Table [5] include the total number of impacts and therefore 
must be considered as upper limits to the glitch events, which will be defined during the 
data reduction. The values of the deposited energy corresponds to the average value 
per hit and < n > to the number of glitches per unit frame rate. It must be noticed 
that the word glitch may not be appropriate here as it does not correspond to a voltage 
shift but to a temperature increase. We keep it in this text but advice the reader to 
consider its meaning as a temperature increase. 
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Table 5 Mean number of impacts per frame rate and the mean heat deposited per hit on 
the blue bolometer grid with spacecraft shielding and secondary events, at solar minimum and 
maximum. 

Blue Grid 



GCR Solar min Solar max 





<n> 




(dT dep )(mK) 


<n> 




(dT dep ){mK) 


H 


(3.9±0.5)xl0" 


-l 


14.0±1.8 


(1.6±0.1)xl0" 


-l 


5.5±0.5 


He 


(0.8±0.1)xl0" 


-l 


6.3±0.5 


(4.3±0.7)xl0" 


-2 


2.8±0.6 


Li 


(1.3±0.1)xl0" 


-3 


0.09±0.01 


(4.9±0.8)xl0" 


-4 


0.04±0.01 


C 


(2.2±0.3)xl0" 


-2 


1.6±0.2 


(0.6±0.1)xl0" 


-2 


0.52±0.02 


N 


(0.7±0.1)xl0- 


-2 


0.52±0.05 


(2.7±0.3)xl0" 


-3 


0.21±0.03 


O 


(3.2±0.2)xl0- 


-2 


2.5±0.2 


(l.l±0.1)xl0- 


-2 


0.9±0.1 



Table 6 Same as Table 5 for the red bolometer grid 

Red Grid 



GCR Solar min Solar max 





(n> 




(dT dep )(mK) 


(n) 




(dT dep )(mK) 


H 


(1.4±0.2)xl0 _ 


-l 


2.9±0.5 


(0.6±0.1)xl0" 


-l 


1.3± 0.3 


He 


(3.9±0.6)xl0- 


-2 


1.5±0.2 


(1.6±0.3)xl0" 


-2 


0.6±0.1 


Li 


(5.1±0.4)xl0 _ 


-4 


0.021±0.001 


(1.8±0.1)xl0" 


-4 


0.009±0.001 


C 


(0.7±0.1)xl0 _ 


-2 


0.37±0.02 


(2.5±0.2)xl0" 


■3 


0.13±0.02 


N 


(2.7±0.5)xl0- 


■3 


0.13±0.02 


(1.0±0.2)xl0- 


■3 


0.07±0.05 


O 


(1.2±0.2)xl0- 


-2 


0.56±0.05 


(3.9±0.4)xl0- 


■3 


0.21±0.02 



7.3 Red bolometer 

We remind that the red bolometer is made of 2 chips and that the grid is larger than 
for the blue one (see No other differences with blue bolometer are present. Here 

we present the results (see Tabled of our simulations normalized to the characteristic 
frame rate of the bolometer (cfr. above) in the 'shield' + secondary events configuration. 

The number of hits between the two arrays differs by a factor 2.6 which is roughly 
the ratio between the difference in pixels number (4) and their volume: 4/1.55. 

7.4 Glitch length 

A large number of glitches (hitting events) are observed on one or two pixels. The 

theoretical glitch length, i.e. the number of nearby pixels affected by an event happening 

in one pixel, can be estimated approximately by simple geometrical considerations on 

the solid angle with respect to the bolometer plane (see Table [7}. 

The distribution of incident angles of the particles hitting the detector is isotropic. 

Some particles will not hit the IPW but will go straight onto the pixels within a solid 

angle, Q\. These ones will not affect the nearby pixels and do not contribute to the 

glitch length. Other particles will hit one or more IPW according to the incident angle 

value, i.e. solid angles J?2, ^3, •■• The incident probability is proportional to the solid 

angle viewed by each pixel and containing all the incoming directions of particles hitting 

the IPW. We may therefore compute the ratio of single pixel to 'two pixels glitches' as 
i?i 

Small corrections to this value are also taken into account: effective cross section for a 
critical angle, approximation of the circular pixel instead of a square one. The related 
uncertainty is given as a lower and upper limit on this theoretical estimation and is 
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Table 7 Comparison between theoretical glitches length (see text for detail) and G4 data for 
the blue bolometer (average from 5 simulations of 9441 protons seed each), without 'shield' 
and without secondary events. In first column the length in pixel units of a glitch. 



Length 


Th up % 


Th; OI1) % 


<G4> % 


2 
3 
4 


9.77 
3.88 
1.80 


6.52 
2.47 
1.13 


14.26l- y 9 



given in Table [7] where these values are compared with the results from the average of 
five proton simulations without a shield. 

We also checked the number of pixels affected by one single GCR for both bolometer 
cameras. We computed the percentage of each GCR type producing glitches affecting 
nearby pixels and assessed that only the adjacent pixels are affected 20 and 30% of 
the time. We also checked the heat propagation effect once a particle hits a IPW: if we 
take the values of the IPW volume (see Appendix A) we find that the maximum heat 
transfer is ^g^-^y) — 0.48 ^/fy along the x-direction and 0.02 jy^y in the y-direction 
(where x-y is the bolometer plane). Both effects may not be negligible and deserves 
further investigation. 

No statistically significant dependencies on solar activity are found. 



8 Observation simulations 

To test the effect of glitches we have simulated observations of a 'strong' and of a 
'weak' source with the PACS Scientific Simulator (Andreani et al., in preparation), a 
complete software code able to simulate from the input sky/source with a given flux 
the entire detection chain up to the output recorded signal. 

We ran two simulations, one with a thermal source at 50 K with a spectral energy 
distribution peaking at 55 fim and very weak one at a temperature of 5 K. 

In the first case the pixels looking at the sources receive an average signal of 3.7 K 
per pixel; a H glitch yields, on average, a temperature increase of 24/3700 i.e. of less 
than 1%. In the case of the very faint source, at 5 K, the heating of each pixel due to 
the source is of 500 mK and the heating produced by the glitch amounts to 24/500, i.e 
~ 5% . The signal can be accurately retrieved only with a suitable observing technique, 
allowing more than one pixel (at least a few) to look at the source, so that the signal 
lost because of the glitch is recovered from the nearby uncontaminated pixels looking 
at the same source. 



9 Test beam results and comparison with simulations 

We report in this section the comparison between the Geant4 simulations and irradi- 
ation tests run in May 2005 at the IPN (Institute de Physique Nucleaire) in Orsay, 
Paris. We ran a complete set of simulations to reproduce the instrumental setup as it 
is described in |10llll| . Details about the simulations can be found in [12] and the full 
description of the instrumental setup with specimen, cryostat and the proton beam are 
defered in Appendix C (see also [1011 11] ). 



12 



Table 8 Comparison of IPN experimental results with the Geant4 output. Column 4 reports 
the flux calculated with the Rutherford scattering formula 13 ; column 5 includes the correction 
for the shadow effect (i.e. the detector sees only 30% of the beam); column 6 gives the measured 
flux; column 7 the flux estimated by Geant4. 



Particles 


Flux 


Au 


Rutherford scat. 


corrected 


measured flux 


Gcant4 




(part/s) 


(nm) 


particles/s 


particles/s 


particle sum 


particle sum 


proton 


76.7 


100 


Ri4.3 


fal.3 


^2.25 


1.86 


proton 


127.7 


100 


Ri7.3 


«2.2 


^3.75 


3.20 



One of the blue chip array was hit by a proton and a particles beams. The number of 
impacts detected in the array and that simulated are reported in Table [8] (see also [11]). 
The different setups during proton irradiations are listed in the first three columns; 
column 4 reports the flux computed with the Rutherford's formula, the fifth the same 
value reduced by 70% because of the filling factor discussed in Appendix C0- Column 6 
lists the measured particle flux and the last one that resulting from Geant4 simulations. 

The experimental Particle sum value is slightly higher than that given by G4. The 
difference can be ascribed to either 1) the cut values used in the simulations for the 
mylar and Al layer placed in front of the detector which may be too large or 2) an 
underestimate of the solid angle for the Rutherford diffusion. A decrease of the cut 
values would increase the number of secondary events, whereas a bigger solid angle 
would increase the primary event rate (proton and alpha seeds). At present it is hard 
to decide which are the most appropriate modifications of simulation parameters. 

Temporal behaviour cannot be easily retrieved from the present simulations and 
in this case we exploit the shadowing effect of the beam cone (see Appendix C and 
Figure[5]) to check for temperature variation as a function of the incoming flux. We have 
compared the pixels which were only slightly affected by the incoming beam (maybe 
even not hit by any particle) with those clearly subject to a strong bombardment. 
We report in Figure U the temperature variation along four imaginary cuts performed 
along the pixel matrix. The temperature varies by at most 15% for pixels receiving the 
largest number of hits. Note that the number of hits on pixel columns 1 to 5 differ 
by a factor > 10 from that on columns 12 to 16. This result was also found in the 
measurements and is due to the geometrical effect of the cone. 

These results show that the particles incoming on a pixel increase its temperature 
by up to 15%. The pixel-to-pixel variation yielded by our simulations, due to the 
increasing intensity of the beam, may be understood in terms of temporal behaviour of 
a single pixel under the effect of a cosmic particle. This may result in a slow responsivity 
drift, i.e. in a change of the response of the single pixel with a time scale depending on 
the bolometer time constant. 

10 Conclusions: effects on observations 

— Red bolometers are less affected by GCR glitches than the blue ones. The pixel grid 
of the red bolometer is a factor of 1.569 larger so that its heat capacity (volumetric 
heat capacity times the volume) is correspondingly higher for the red pixels. The 
number of impacts is proportional to the detector surface. Both the red and the 

5 this filling factor was introduced to take into account the experimental geometry, the 
bolometer itself and the grid. 
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Fig. 4 Temperature variation per unit time (Is) under proton irradiation along four sections 
of the 16 X 16 pixel matrix, after subtracting the blind pixels: the two median sections (vertical: 
x pixel =10, upper panel; horizontal y v i xe i=&, second panel from the top) and the two diagonal 
sections (down-up: y p i X el =x pixer^j third panel from the top; up-down: y p i xe i=18-x p i xe i bot- 
tom panel). In black the data referring to the experiment 3 setup, in red for exp. 4, in green 
for exp. 5 and in blue for exp. 6. The shadow effect is less present in the plot referring to the 
up-down diagonal, whose pixels suffer almost the same shadowing from the cone. 
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blue bolometer are affected by the same impinging flux but they have different 
collecting areas: 2048 grids for the blue and 512 grids for the red and the expected 

ratio of hit rate is: ^"'J^ 8 = ^2392 ^ ^ltf ~ 2 ' 55 ' where ^bluo and A rcd are 
the array areas of the blue and red bolometers respectively. 

— We have quantified the contribution of secondary events produced by the 'shield'. 
This latter acts as the major source of secondary events, while the bolometer itself 
produces approximately 1/3 of them. 

— The sensitivity of the bolometer cameras may be limited by glitch impacts for very 
faint sources; we expect a glitch contamination < 5-6% of the background. The 
contamination is lower for the red camera. 

— A test with the instrument simulator shows that in case of very faint sources a 
suitable observing technique and glitch removal algorithm must be put in place. In 
particular this latter should work on each chopping time interval in order to avoid 
false detections due to the sudden change of the flux level. 
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11 Appendix A: The Detector Design 

In the G4 code a class describing the geometry and chemical composition of the detec- 
tors must be set. Here we describe how the geometrical model for the PACS bolometers 
is built. 

The detector module is made of different parts: a grid, an Inter-Pixel- Wall, differ- 
ent layers (Al, Si, Ti), a Si strip, reflecting slabs, etc. [3J. Each pixel is supported by 
a silicon grid 5fim thick with a thermometer, and is suspended in an Inter-Pixel- Wall 
(IPW) grid 450/im thick which has a reference thermometer and acts as a thermostat. 
The detector is 'decomposed' in different cells, each of a simple geometrical shape (par- 
allelepipeds (called boxes by G4) and spheres). Each part is characterized by its shape, 
its chemical composition and spatial position (in a Cartesian 3D space coordinates). 
The output signal is a function of the temperature difference between the pixels and 
the reference thermometers. 

Attached to them two columns of blind pixels heated at a mean temperature of nom- 
inal pixels are used to suppress any undesired electronic temporal deviations. All the 
components of the detector are at a temperature of 0.3 K. 

— The four aluminum layers have been approximated as boxes 0. One large layer 
(43x50.88x2.5 mm) placed at the central spatial coordinate system (0,0,0). A 
second small Al layer (37.95x50.88x2.5 mm) lying horizontally on the large layer. 
Two more layers, vertically placed and tangentially to the x base of the two previous 
layers (43x1.56x9.7 mm). 

— The bolometer camera is made of two different arrays: the blue one consists of 4x2 
chips and the red one of 2x1 chips. Each chip contains the following components 
(single or multiple): 

1. One Titanium Layer : an L-shaped volume which is simulated with two boxes 
lying one under the other. The large box (21.7x12.2x2 mm) is placed 0.4 mm 
above the small Al layer. Every Ti box is detached by the other by 0.55 mm. 
The smaller Ti box (we called it 'Ti edge', 3x12.2x1.9 mm) is lying below and 
on the external parts of the large Ti layer. 

2. The Silicon Layer (Si Layer) is a unique box (19.15x12.63x0.45 mm) placed 
0.8 mm above the large Ti layer. Each Si layer is detached one by the other by 
0.12 mm. 

3. the 16x18 pixels of each chip have three identical Reflecting Slabs made of 
Titanium, Palladium and Gold. They lie one on top of each other over the Si 
layer and have dimensions of 640x640x0.3 /im 

4. Each pixel has a Silicon grid which is approximated by a continuum slab with 
a volume equivalent to the real one and a thickness equal to the real one. It 
turns out that the slab has a projected (x,y) area smaller that the effective 
grid. We draw it as a box of dimensions 416.653x416.653x5 /jm, for the blue 
bolometer, and (521. 912x521. 912x5/^m), for the red bolometer. The volumes 
of the absorbers are then 868000 /im 3 for the blue and 1361960 /im 3 for the 
red. 

Analysis of the different covering factor of the slab, compared to the real grid 
will be given below. 

5. 17x 19 Indium spheres (radius 10 pm) are lying on top of the Si layer. They are 
placed at the intersection of the Inter-Pixel- Wall (IPW), sustaining it. 



In case of a box shape we gave its measures as xxyxz. 
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6. The Inter-Pixel- Wall made of Si defines the border of the Si grid. The IPW, 
pole-on viewed, is a grid and it was decomposed in similar parts: 5 different 
boxes are built making the whole IPW. 2+2 boxes make the border of the 
IPW thicker respect to the internal IPW (the ones separating the grids): (1) 
the longest boxes (horizontal border, HBWall) are 14.13x0.35x0.45 mm; (2) 
the shortest boxes (vertical border, VBWall) are 0.35x11.93x0.45 mm. (3) The 
external vertical border (Extra Wall) is thicker than the usual VBWall. That 
is we added an extra box of 1.65x12.13x0.45 mm. 

The IPW grid was designed by considering (4) 15 long (horizontal, HWall) 
boxes of 13.41x0.07x0.45 mm each, with (5) 17 short boxes (vertical, VWall) 
placed vertically to the long ones (0.07x0.68x0.45 mm each). 

7. The Mask (blind pixels, 2x16). 

Both the bolometers have two rows of blind pixels. They are equal to the other 
pixels, but front masked by different strata of different materials. These blind 
pixels are placed in the external part of the chip. There are three slabs of the 
same material and in the same order of the reflecting layers, that is Ti, Pd and 
Au, designed as boxes 3.125mmx 12.63mm x0.3/im each. Above them there is 
a thick Silicon 3.125x12.63x0.45 mm layer. 



12 Appendix B: The Physics Model 

We define the list of physical interactions with the shielding and the different compo- 
nents of the photoconductors in order to best reproduce the impact of GCRs on the 
photoconductors. The list is comprehensive of all the possible physical interactions: 

— Electro-magnetic physics: photo electric and Compton effect; pair production; elec- 
tron and positron multiple scattering, ionization, bremsstrahlung and synchrotron; 
positron annihilation. 

— General physics: decay processes. 

— Hadron physics: neutron and proton elastic, fission, capture and inelastic processes; 
photon, electron and positron nuclear processes; 7r + , 7r~, K + , K — , proton, anti- 
proton, E + , E~ , anti-X^, anti-X^, E + , 2~ , anti-S , anti-£", Q~ , anti- 
Q + and anti-i?~ multiple scattering and ionization. We added the High Precision 
Neutron dataset, that is valid for neutrons till 19.9 MeV. 

— Ion physics: multiple scattering, elastic process, ionization and low energy inelas- 
tic processes. We added inelastic scattering (J. P. Wellish, private communication) 
using the binary light ion reaction model and the Shen cross section [14] , 

— Muon physics: fi + , fx~ multiple scattering, ionization, bremsstrahlung and pair 
production; pT capture at rest; r + and r _ multiple scattering and ionization. 

For 'ion physics' of aZp/ia-particles, protons, and H isotope nuclei we used the cross 
section by Tripathi et al. ([H]), holding for energies of up to 20 GeV per nucleon. The 
cross sections used for the other ions |16| are valid only up to 10 GeV per nucleon. 
This means that we had to cut out high energy ions of the GCRs (whose flux is lower 
than 10~ 5 /m 2 /s/sr). The treatment of ion physics is strictly valid only for nuclei not 
heavier than C. However we have applied it also to simulate N and O, which should 
not behave much differently. 
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The model called by our physics list is different depending on the energy of the 
incoming particles. Our physics list is based on QGSP_BIcfl list (J.P. Wellish, private 
communication; see the Geant4 web page User Guide §2.4.3 for details). Modifications 
to such a list are: 

1. We added the High Precision Neutron dataset, which is valid for neutrons up to 
19.9 MeV. 

2. For neutrons and protons, we introduced^ the Bertini cascade [T7] up to 90 MeV 
and made the BIC starting from 80 MeV. 

3. QGSP is assumed valid from 10, rather than 12, GeV. 

4. We have taken into account the gamma- and electro-nuclear reactions, so we con- 
sidered the electromagnetic_GN physics list. 



13 Appendix C: the description of the set up of the irradiation 
measurements 

The detector used for the proposed investigations under proton and alpha particle ir- 
radiation was the detector module SMD 3 blue (N° 05-05 type MV), i.e. the Flight 
Model (equivalent doping) mounted on a common carrier in flight configuration. The 
"3" indicates the position on the common carrier. This means that the mounting is as 
if the bolometer was complete, so 7 locations out of 8 on the mount were empty. Both 
the bolometer and the reaction chamber were placed in vacuum. Two pumps held the 
vacuum to ~ 10~ 6 mbar. 

The proton energy in the experimental setup was 20 MeV and the flux could be changed 
from 76.6 to 127.7 particles per second. 

The beam, before reaching the blue bolometer elements penetrates several layers and 
undergoes Rutherford diffusion due to a 45° inclined gold foil placed at lm from the de- 
tector, decreasing the flux of the beam. The foil has two different thicknesses according 
to the running test. Proton tests were all run with a 100 nm thick foil. 

The diffused beam then encounters a couple of 14/im thick Al foils, placed 75 cm 
from the detector. Then it crosses a 45° inclined foil of mylar (a polymer: [-0-C::0- 
C6H6-C::0-0-CH2-CH2-]n) 25 /^m thick; it goes through a sort of cave cone. Then 
the beam hits the blue bolometer module. Because of the detector geometry - the 
bolometer itself does not consist of a full Si chip but of a grid - the surface exposed to 
the cone is not 100 percent full. In addition the cone itself introduces a shadow effect 
on the bolometer because the cone face-on area masks part of it. This can be seen 
in Figure [5] where the number of glitches per second and pixel are shown in a three 
dimensional plot. A shadow effect is present on the left, the bolometer array receives 
less hits because of the cone. Considering the surface of the bolometer array only 30% 
"sees" the beam. 



7 Quark Gluon String Physics_Binary Cascade 

8 It is well known that (1) the BIC does not reproduce well data for energies below 50 
MeV I16| . In contrast, the Bertini cascade works well below 50 MeV for all but the lightest 
nuclei. And (2) the BIC for pions and kaons has not been tested yet (J.P. Wellish, private 
communication) . 
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Fig. 5 The picture shows the number of observed glitchcs/s/pixcl. On the left the shadow 
due to the cone is visible. Considering the surface of the bolometer arrays (16*16 pixels or 
bolometers) only 30% of it is "filling" (due to geometrical configuration, see text for details). 



